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galaxies and implications for small-scale fluctuations at 
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ABSTRACT 

This paper investigates the clustering properties of SCUBA-selected galaxies in the 
framework of a unifying scheme relating the formation of QSOs and spheroids (Granato 
et al. 2000). The theoretical angular correlation function is derived for different bias 
functions, corresponding to different values of the ratio iV/haioZ-^sph between the mass 
of the dark halo and the final mass in stars. SCUBA sources are predicted to be 
strongly clustered, with a clustering strength increasing with mass. We show that 
the model accounts for the clustering of Lyman-break Galaxies, seen as the optical 
counterpart of low- to intermediate-mass primeval spheroidal galaxies and is also con- 
sistent with the observed angular correlation function of Extremely Red Objects. Best 
agreement is obtained for Mhaio/Afsph = 100. We also consider the implications for 
small scale fluctuations observed at sub-mm wavelengths by current or forthcoming 
experiments aimed at mapping the Cosmic Microwave Background (CMB). The pre- 
dicted amplitude of the clustering signal in the 350 GHz channel of the Planck mission 
strongly depends on the halo-to-bulge mass ratio and may be of comparable amplitude 
to primary CMB anisotropics for multipole numbers I > 50. 

Key words: galaxies: clustering - galaxies: infrared - cosmology: theory - large-scale 
structure - cosmology: cosmic microwave background 



1 INTRODUCTION 

The advent of the 450-850 ^ra Submillimetre Common User 
Bolometric Array (SCUBA) camera (Holland et al. 1999) 
has opened the possibility of exploiting the enormous poten- 
tial of sub-mm astronomy for investigating the early phases 
of galaxy evolution, taking advantage of the strongly nega- 
tive K-correction which boosts the visibility of high-z galax- 
ies (Smail, Ivison & Blain 1997; Barger et al. 1998; Hughes 
et al. 1998; Barger, Cowie & Sanders 1999; Bales et al. 1999; 
Blain et al. 1999a). 

A variety of models has been proposed to reproduce 
the properties of the submillimetre population (e.g. Guider- 
doni et al. 1998; Blain et al. 1999b; Devriendt, Guiderdoni 
& Sadat 1999; Devriendt & Guiderdoni 2000; Takeuchi et 
al. 2000; Pearson 2000), in particular the dramatic star- 
formation rate inferred from the data. 

Granato et al. (2000 - hereafter Paper I) propose a uni- 
fied scheme for the formation of the QSO and spheroidal 
population. According to their model, at the time of the on- 



set of the QSO activity, the host galaxy has already under- 
gone a major star-formation event, with a characteristic time 
Tb ranging from ^ 0.5 to ~ 2 Gyr when going from more 
to less massive hosts. The QSO onset practically stops the 
star-formation process, and when eventually the reservoir of 
gas fueling the Active Galactic Nuclei is exhausted, the host 
galaxy will evolve in a passive way, with a spectrum becom- 
ing quickly red. This approach identifies SCUBA galaxies as 
the progenitors of QSO hosts and provides an excellent de- 
scription of (amongst others) the observed 850 /xm number 
counts. 

In this paper we present predictions on clustering prop- 
erties of SCUBA galaxies in the framework of the Granato 
et al. (2000) model. Such properties are determined by the 
combination of the evolution of the underlying mass fluc- 
tuations and the bias relating galaxy overdensities to mass. 
Since the bias factor is determined by properties related to 
the process of galaxy formation, it follows that measure- 
ments of clustering can lead to very precious insights on the 
galaxy population producing the signal. 
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Unfortunately the number of SCUBA galaxies observed 
is still too small to allow for any clustering measurement. 
The model by Granato et al. (2000) however establishes a 
tight link between this population and Lyman-break galax- 
ies (Steidel et al. 1996; LBGs hereafter), interpreted as 
lower-mass/lower-submm luminosity sources. Such model 
can thus be tested against measurements of LEG cluster- 
ing (Giavalisco et al. 1998). A further test is provided by 
the recent detection (Daddi et al. 2000) of strong cluster- 
ing of Extremely Red Objects (EROs), which were found to 
be the optical counterparts of two of the still few SCUBA 
sources having reliable optical identifications (Small et al. 
1999). 

An important issue related to the above discussion is the 
effect of clustering on small scale anisotropies at mm/sub- 
mm wavelengths, where experiments to map the Cosmic Mi- 
crowave Background either have been (e.g. BOOMERanG, 
MAXIMA, TopHat, ACBAR, ARCHEOPS) or are about 
to be carried out (notably the ESA's Planck mission). Small 
scale fluctuations of the background intensity due to cluster- 
ing of unresolved sources are proportional to the amplitude 
of their two-point correlation function and to their volume 
emissivity (see, e.g., De Zotti et al. 1996). The deep SCUBA 
surveys have demonstrated that the bulk of the extragalactic 
background at 850 fim comes from a relatively small number 
of sources (at least ten time less than those accounting for 
the bulk of the optical/UV background) at substantial red- 
shifts. Combining the (rather limited) direct spectroscopic 
information (e.g. Prayer 2000) with estimates based on the 
sub-mm to radio spectral index (Carilli & Yun 1999), it is 
concluded that the median redshift of the sub-mm popula- 
tion brighter than 1 mjy is probably ~ 2.5-3 (Small et 
al. 2000). It then follows that most of the background comes 
from a population of extremely luminous (hence presumably 
very massive) high- 2: galaxies which probably correspond to 
the rare, high density peaks in the primordial matter distri- 
bution. The standard structure formation scenario predicts 
that these peaks are strongly clustered (Kaiser 1984; Baron 
& White 1987). Also, according to the currently favoured 
view, SCUBA sources probably correspond to the phase of 
intense star formation in large spheroidal galaxies, a process 
that appears to be completed for 2; > 1; as a consequence, 
the relevant redshift range is limited and dilution of the clus- 
tering signal is relatively small. 

The amplitude of intensity fluctuations due to cluster- 
ing is correspondingly expected to be large. Scott & White 
(1999) worked out a preliminary estimate adopting for sub- 
mm sources the angular correlation function determined for 
z ^ 3 LBGs (which however are mostly undetectable by 
SCUBA; see Chapman et al. 2000). Here we provide quan- 
titative estimates of this effect using a technique developed 
by Peacock & Dodds (1996), Matarrese et al. (1997) and 
Moscardini et al. (1998), which takes into account the ef- 
fects of the past light-cone, the linear and non-linear growth 
of clustering and the evolution of the bias factor. The epoch- 
dependent mass distribution of sources is modelled following 
Granato et al. (2000). 

The plan of the paper is as follows. Section 2 will be 
devoted to the analysis of the clustering of SCUBA galax- 
ies and comparisons with the observed clustering properties 
of LBGs and EROs, while in Section 3 we use the findings 
of Section 2 to make predictions on the temperature fiuctu- 



ations produced by unresolved SCUBA-like sources and to 
discuss their implications for experiments aimed at mapping 
the CMB at sub-mm wavelengths. Section 4 summarizes our 
conclusions. 

Throughout the paper we will assume a cosmological 
model with ho = 0.7, Slo = 0.3, A = 0.7, and adopt the nor- 
malization of the rms mass fiuctuation as ~ 1 determined 
by the 4-year COBE data (Bunn & White 1997). 



2 CLUSTERING OF SCUBA SOURCES 

The angular two-point correlation function ■w{6) is defined as 
the excess probability, SP, with respect to a random Poisson 
distribution, of finding two sources in the solid angles SQi 
S^l2 separated by an angle 9, and it is defined by 



SP = nsniSn2 [1 + w{e)] 



(1) 



where n is the mean number density of objects in the par- 
ticular catalogue under consideration. 

The theoretical expression for w{0) can be derived from 
its spatial counterpart ^ by projection via the relativistic 
Limber equation (Peebles 1980): 

_ /o°°/o°° N^z) bgff (M„,in, z) (dz/dx) gr, z) dz du 



dz\ 



where x is the comoving radial coordinate, r — {u^+x^O^y^'^ 
(for a flat universe and in the small angle approximation: 
6 <C Irad), and A^(2:) is the number of objects within the 
shell (2, z + dz). 

The mass-mass correlation function S,{r,z) to be in- 
serted in eq.(^) has been obtained following the work by 
Peacock & Dodds (1996, but also see Magliocchetti et al., 
2000), extended by Matarrese et al. (1997) and Moscardini 
et al. 1998) to redshifts z > 0, which provides an analyti- 
cal way to derive the trend of ^(r, z) both in the linear and 
non-linear regime. Note that ^(r, z) only depends on the 
underlying cosmology and on the normalization of erg. The 
relevant properties of SCUBA galaxies are included in the 
redshift distribution of sources N{z), where N{z) dz gives 
the number of objects in the shell (2, 2-|-d2) (given in Paper 
I and shown in Fig. ^ for fluxes respectively greater than 1, 
10, 50 mJy), and in the bias factor &cff (Afmin, z). 

The effective bias factor 6cff(Mmin, 2) of all haloes with 
masses greater than some threshold Mmin is obtained by 
integrating the quantity b{M, z) - representing the bias of 
individual haloes of mass M - opportunely weighted by the 
number density nscuBA{M, z) of SCUBA sources: 



bcS (2) 



'^^ ^(^^' '*scuba(M, 2) 
dM nscvBA{M,z) 



(3) 



Note that as nscuBA can be thought as the fraction of haloes 
hosting a galaxy in the process of forming stars, its expres- 
sion can also be written as nscuBA^M, 2) = n(M, 2) Ts/th, 
where n{M, 2) dM is the mass spectrum of haloes with 
masses between M and M + dM (Press & Schechter, 1974; 
Sheth & Tormen, 1999), Tb is the duration of the star- 
formation burst and th is the life-time of the haloes in which 
these objects reside (see Martini & Weinberg, 2000). The du- 
ration of the star-formation phase Tb varies with the mass 
of the spheroid undergoing the burst (see Paperl), assuming 
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Figure 1. Predicted redshift distribution for SCUBA galaxies re- 
spectively brighter than 1, 10 and 50 mjy. The Figure is adapted 
from Paper I. 



a value Tj, — 0.5 Gyr in the case of massive objects (i.e. with 
bulge masses Msph ^ 1-5 IO^^Mq) and increasing when go- 
ing to lower masses. 

We take the expression for b{M, z) to be plugged into eq.(^ 
from Jing (1998): 



6(M,z)= ( 1 + - 



a2(M) 



a-^(M) 
2 5i{z) 



+ 1 



(6-^''cff) 
100 



(4) 



with 5c{z) — 5c/D{z), D{z) being the linear growing factor 
and Sc the critical density for collapse of a homogeneous 
spherical perturbation {5c = 1.686 in an Einstein-de Sitter 
universe). a^{M) is the variance of density fluctuations on 
the mass scale M, and rieff = — 3 — 6(dln(T/cilnM). Note that 
the assumption of a linear bias is justified by the redshift 
distribution of the sources, all appearing for z > 0.75 (see 
e.g. Somerville et al., 2001). 

According to Granato et al. (2000), sources showing up 
in the SCUBA counts can be broadly divided into three cate- 
gories: low-mass (masses in the range M^ph — — 10^" Mq, 
duration of the star formation burst Tb ~ 2 Gyr, and typ- 
ical fluxes S < 1 mJy), intermediate-mass (Msph — 10^"- 
1O"M0 and Tb ~ 1 Gyr) and high-mass (M^ph > IO^Mq, 
Tb ~ 0.5 Gyr, dominating the counts at fluxes S > 5 — 
10 mJy). Note that by Msph we denote the mass in stars at 
completion of the star formation process. 
In order to evaluate the bias factor in eq.(^ we then consider 
two extreme cases for the ratio between the mass in stars 
and the mass of the host dark halo: Mhaio/Msph = 100 and 
Mhaio/Msph = 10. Note that Mhaio/M^ph = 10 roughly cor- 
responds to the ratio Slo/f^baryon between total and baryon 
density, where we adopted for the latter quantity the stan- 
dard value from primordial nucleosynthesis; this corresponds 
to having assumed all the baryons to be locked into stars 
and, as a consequence, gives a conservative lower limit to 
the quantity Mhaio/Mgph. Mhaio/Mgph — 100 is instead re- 
lated to ilo /^* , being the present mass density in visible 



Figure 2. Predictions for the angular correlation function w{9) 
for different halo-to-bulge mass ratios and different flux cuts. Solid 
lines are for sources brighter than 50 mJy, dashed lines for sources 
brighter than 10 mJy, while dotted lines correspond to 5 > 1 mJy. 
Higher curves of each type correspond to Mhaio/^^sph = i^OO, 
lower ones to A/jjaio/^'^sph = 10. 



stars (Persic & Salucci, 1992; Fukugita et al. 1998). The 
likely value should be Msph/Mhaio = 1 — 3 % (cf. Paper I). 

Armed with the above results we can finally evaluate 
the two-point correlation function in eq. (^) for different 
Afhaio/Afsph ratios and different flux cuts. Figure ^presents 
our predictions for w{8), respectively for a flux cut of 50 
(solid line) , 10 (dashed line) and 1 (dotted line) mJy. Higher 
curves of each kind correspond to the case Mhaio/Msph = 
100, while lower curves refer to the Mhaio/Afsph = 10 one. 

The highest clustering amplitude is found for the bright- 
est sources (S > 50 mJy). This is because they are associ- 
ated to the most massive dark halos and are therefore highly 
biased tracers of the dark matter distribution. In addition, 
according to the model by Granato et al. (2000), they have 
a rather narrow redshift distribution (see Fig. |l|) , so that the 
dilution of the clustering signal is minimum. 
The very sharp drop of all the curves at ~ 1° is due to the 
absence of nearby (z < 1, see Fig. |l|) objects. This reflects 
the notion that the actively star-forming phase in spheroids 
is completed at z > 1. 

The model introduced in Paperl establishes a strong 
link between SCUBA sources and LBGs where this latter 
population is predicted to be the low- to intermediate-mass 
tail of primeval spheroidal galaxies, with Tb ~ 1 — 2 Gyr 
and a star-formation rate ranging from a few to a hundred 
Mq yr-i. 

Measurements of the angular correlation function of LBGs 
(Giavalisco et al. 1998) then offer a way to test such model. 
In Fig. ^ we plotted the expected w{6) for sources with 
S > 1 mJy (corresponding to Msph ^ 10^" Mq) within the 
redshift range 2.5 < z < 3.5, covered by the Giavalisco 
et al. (1998) sample. As in Fig. the higher curve is for 
Mhaio/Msph = 100, the lower one for Mhaio/Afsph = 10. The 
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Figure 3. Angular correlation function of LBGs. Model pre- 
dictions are shown by the dotted lines (the higher one is for 
Mjjaio/^sph = 100, the lower for A/jjaio/A^sph = 10), while data 
points represent the Giavalisco et al. (1998) measurements. 



data points show the Giavalisco et al. (1998) measurements. 
Despite the large errors, it is clear that the predicted trend 
for w{6) correctly reproduces the data, a better fit obtained 
for Mhaio/Msph = 100. This result, apart from providing a 
further proof of the tight link between SCUBA galaxies and 
LBGs, also confirms the expectations of Paper I for a small 
fraction (on the order of a few percent) of the total mass to 
be confined into stars. 

Strong clustering of EROs has been recently detected by 
Daddi et al. (2000). This result may again provide a relevant 
test for our model since it seems likely for a substantial frac- 
tion of SCUBA galaxies to be EROs with R~K > 6: at least 
two out of nine (> 20%) sources in the Cluster Lens Survey 
belong to this class (Small et al. 1999), while other submil- 
limeter surveys have discovered additional EROs as coun- 
terparts of SCUBA sources (Ivison et al. 2000). For their 
ERO sample with Ks < 19.2 and R - Ks > 5.3, Daddi et 
al. (2000) find w(6l) ~ 0.07 for 6 ~ 0°.l, not far from the 
predictions in Fig. 



3 SMALL SCALE FLUCTUATIONS FROM 
DISCRETE SOURCES 

3.1 The Intensity Correlation Function 

An issue intimately connected with the analysis of galaxy 
clustering performed in the previous Section is the study 
of the contribution of unresolved sources (i.e. sources with 
fluxes fainter than some detection limit Sd) to the back- 
ground intensity. Its general expression is: 



dN 

Is 



S dS 



= -y^.logLL 



dz $(L, z] 



K{L,z) dx 

{1 + zy dz ^ ' 



(see e.g. De Zotti et al. 1996), where dN/dS denotes the 
differential number counts, Lmax and Lmin are respectively 
the maximum and minimum local luminosity of the sources, 
K{L, z) — {l + z)L[i'{l+z)]/L{i') is the K-correction, Zmax is 
the redshift when the sources begin to shine, z{Sd, L) is the 
redshift at which a source of luminosity L is seen with a flux 
equal to the detection limit Sd, ^{L,z) is the luminosity 
function (i.e. the comoving number density of sources per 
unit dlogL), and x is the comoving radial coordinate. 

The intensity fluctuation 51 due to inhomogeneities in 
the space distribution of unresolved sources is then given by 
eq. (^, with the quantity $(!/, z) replaced by 5${L, z). It is 
easily shown that the angular correlation of such intensity 
fluctuations 



c{e) = {5i{9',^') 5i{e"A")), 



(6) 



where {d',4'') and [9" ,(f)") define two positions on the sky 
separated by an angle 6, can be expressed as the sum of two 
terms Cp and Cc, the first one due to Poisson noise (i.e. 
fluctuations given by randomly distributed objects), and the 
second one owing to source clustering. 

Scott & White (1999) have shown that if dusty galax- 
ies cluster like LBGs (Giavalisco et al. 1998), anisotropies 
due to clustering dominate the Poisson ones at all angular 
scales. According to our model, LBGs are the optical coun- 
terparts of the low- to intermediate- mass tail of SCUBA 
galaxies; as illustrated in Figure ^, this makes the clustering 
signal of LBGs a factor ~ 10 lower than the one obtained for 
SCUBA galaxies; therefore we can safely assume the Poisson 
term Cp to be negligible in the following analysis and only 
concentrate on temperature fluctuations caused by the Cc 
term (hereafter simply called C). 

By making use of the quantities defined in Section 2 
and of eq. (^, the clustering term in eq. (^) takes the form: 



{1 + zY 



- duiir 



(7) 



where 

" ${L,z) K{L,z) LdlogL (8) 

is the effective volume emissivity (see Toffolatti et al. 1998). 

Following the prescriptions of Section 2, we have then 
evaluated C{6) in eq.(^) separately for the three cases of 
low-, intermediate- and high-mass objects, by plugging in 
eq. (pi) the appropriate expressions for the luminosity func- 
tion. Note that in this case, the quantity 6^(Mmin, z) should 
indeed be read as b^(Mmin, Mmax, 2), where Mmax is the 
maximum halo mass corresponding to the maximum visi- 
ble bulge mass (i.e. upper limit for the mass locked into 
stars whose values are given in Section 2) . This corresponds 
in eq. (^) to a replacement in the upper limit of the integrals 

of 00 with Afmax. 

Figure ^ shows predictions for the correlation signal due 
to intensity fluctuations (in units ergs^ cm"** Hz~^ sr~^) 
in the case of SCUBA galaxies. Dotted, dashed and solid 
lines are respectively for spheroids of M™Jj" = 10^, 10^" 
and 10^^ Mq. In each plot the two solid curves correspond 
to different detection limits, the 100 mjy case represented 
by the higher curve and the 10 mJy case by the lower one. 
Note that, due to the steepness of the bright end of source 
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Figure 4. Left panel: predictions for the correlation signal due to intensity fluctuations of SCUBA galaxies in the case Mhaio/A^sph = 10. 
Dotted, dashed and solid lines are respectively for spheroids of low (M^ph = 10^-10^" Mq), intermediate (Afsph = lO^^-lO^^ ^q), and 
high (Mgpij > 10^^ ^q) masses. The two solid curves are obtained for different flux detection limits, the 100 mjy case represented by 
the higher curve and the 10 mJy case by the lower one. Right panel: same as before but with a ratio Af ij^io /A^sph = 100. 



counts, the results are independent of Sd if Sd > 100 mJy. 
Results are shown for Mhaio/A/gph = 10 (left-hand panel) 
and Mhaio/Afsph = 100 (right-hand panel). Note once again 
the steep drop of all the curves around ~ 1° caused by the 
absence of nearby sources of this kind and the very small 
contribution given by the low-mass objects, mostly due to 
their small bias function and their modest contribution to 
the effective volume emissivity. 

In order to evaluate the total contribution of the clus- 
tering to intensity fluctuations one has to add up all the 
values of C{9) in eq. (^) obtained for the different mass in- 
tervals and also to take into account the cross-correlation 
terms between objects of different masses, so that the final 
expression for C{9) is given by 



(9) 



where the indexes i,j stand for high, intermediate and low 
masses. 



3.2 Power Spectrum of Temperature Fluctuations 

The angular power spectrum of the intensity fluctuations 
[eqs. (^-(^)] can be obtained starting from the well known 
expression for the expansion of temperature fluctuations in 
spherical harmonics 



ST 



(10) 



If the fluctuations are a stationary process it can be shown 
(Peebles 1993) that the previous expression does not depend 



on the index m (i.e. on the 4> angle), so that we can write: 



C, 



21 + 1 ^ 

m — — l 



with 





0.1 = 



T 



ST{9) 21 + 1 



47r 



PiicoaO) sine (16 ( 



(11) 



(12) 



'0 Jo 

(see Toffolatti et al. 1998), where Pi{cos9) are the Legendre 
polynomials, and ST/T is expressed as 



ST 



,WH(f)- 



^1/2 ^ 



AYC™T(6I) 



2 ktT 



exp 



KhTj 



X exp 



(- 



kbT 



(13) 



which relates intensity and temperature fluctuations 
(T=2.726 K Mather et al. 1994). 

Fig ure |5| shows t he predicted values for the quantity 
STi = + l)Ci/2n (in units of K) at 353 GHz (850 fim) 
- the central frequency of one of the channels of the High 
Frequency Instrument (HFI) of the ESA's Planck mission - 
as a function of the multipole I up to I — 1000. Results are 
plotted for two values of the source detection limit {Sd = 100 
and 10 mJy; note that the source detection limit for this 
Planck channel is likely to be > 100 mJy) and the usual 
two values of Mhaio/Aisph. Also shown, for comparison, is 
the power spectrum of primary (CMB) anisotropies (solid 
line) predicted by a standard Cold Dark Matter model for 
a ACDM (A = 0.7, Qo = 0.3, ho = 0.7), computed with the 
CMBFAST code developed by Seljak & Zaldarriaga (1996). 

At this frequency our model predicts fluctuations of am- 
plitude due to clustering comparable to (and possibly even 
larger than) those obtained for primary CMB anisotropies 
at I > 50. As illustrated by Fig. M, most of the clustering 
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Figure 5. Predicted power spectrum of temperature fluctuations 
STi = ^Jl{l + 1)Ci/2tt (in units of K) as a function of the multi- 
pole / at 353 GHz, the central frequency of a Planck/HFI channel. 
The dashed lines are for a detection limit Sd = 100 mjy, the dot- 
ted ones for Sd = 10 mJy. In both cases the higher curves are ob- 
tained for Miiaio/^4ph=100, the lower ones for Miiaio/Afsph=10. 
The solid line represents the power spectrum of primary CMB 
anisotropies as predicted by a standard Cold Dark Matter model 
for a ACDM cosmology (A = 0.7, = 0.3, ho = 0.7). 

signal comes from massive galaxies with fluxes S > 1 mJy, 
which lie at substantial redshifts (see Fig. ^ and are there- 
fore highly biased tracers of the underlying mass distribu- 
tion. Also the strongly negative K-correction increases their 
contribution to the effective volume emissivity [eq. (^)] and 
therefore to the fluctuations. 

We therefore expect this effect to give a substantial con- 
tribution to the power spectrum extracted e.g. from the 
BOOMERanG map (de Bernardis et al. 2000) at 400 GHz. 
For I > 200 its contribution may be comparable to that due 
to inhomogeneous emission by the interstellar dust. 



4 CONCLUSIONS 

We have worked out predictions for the clustering properties 
of SCUBA-selected galaxies in the framework of an astro- 
physically grounded model, relating the formation of QSOs 
and spheroids (ellipticals, SO galaxies and bulges of spirals). 
The theoretical angular correlation function has been de- 
rived for different bias functions, corresponding to differ- 
ent values of the ratio Mhaio/Msph between the mass of a 
spheroid locked in stars and the mass of its host halo. 

SCUBA sources are predicted to be strongly clustered, 
with a clustering strength which increases with increasing 
mass (or equivalently with increasing luminosity), as a con- 
sequence of the fact that they are very massive and shin- 
ing at substantial redshifts. Since the clustering amplitude 
strongly depends on the quantity Mhaio/A^aph, future mea- 
surements of the angular correlation function w{9) will be 
able to discriminate amongst different models of SCUBA 



galaxies and in particular to determine the amount of bary- 
onic mass actively partaking the process of star formation. 

Under the hypothesis of low- to intermediate-mass 
primeval spheroids to show up in the optical band as LBGs 
(as argued by Granato et al. 2000), we have then compared 
our predictions with the clustering measurements obtained 
by Giavalisco et al. (1998) for a wide sample of LBGs at 
z ~ 3. The agreement is good for a mass ratio between the 
dark halo and the stellar component in the range 10-100; 
high values for this ratio seem to be favoured. 
Also, the predicted amplitude of the angular correlation 
function of SCUBA galaxies is consistent with the one de- 
termined by Daddi et al. (2000) for EROs with R-Ks> 5.3 
and Ks < 19.2. This would support the observational evi- 
dence for a substantial fraction of SCUBA galaxies to be 
identified with EROs (Small et al. 1999; Ivison et al. 2000). 

We have also considered the effect of clustering on maps 
of the Cosmic Microwave Background at sub-mm wave- 
lengths, with special attention to the 353 GHz (850 /xm) 
channel of the High Frequency Instrument (HFI) of the 
ESA's Planck mission. The possibility of large ffuctuations 
due to clustering at this (and also shorter) wavelength was 
first pointed out by Scott & White (1999) who adopted 
a phenomenological approach, assuming for SCUBA galax- 
ies the clustering amplitudes observed for LBGs (Giavalisco 
et al. 1996). Our treatment, leading to a derivation of 
the angular autocorrelation function in the framework of 
the currently standard hierarchical clustering scenario, indi- 
cates that massive SCUBA galaxies (expected to dominate 
the population of sources sampled at fluxes S > 10 mJy) 
are likely to be even more strongly clustered than LBGs. 
The corresponding fluctuations are then expected to be, at 
~ 350 GHz, of an amplitude comparable with those due to 
primary CMB anisotropies for multipole numbers I > 50. 
This would imply that important information on the cluster- 
ing properties of faint sub-mm galaxies (and hence on their 
physical properties such as their mass and/or the amount 
of baryons involved in the star-formation process) resides in 
the BOOMERanG maps at 400 GHz where, however, the 
dominant signal is expected to come from interstellar dust 
emission. On the other hand, if the ratio Mhaio/Afsph is large, 
we expect the clustering signal to be comparable with dust 
emission fluctuations at I > 200. It might therefore be possi- 
ble to separate these two contributions and possibly exploit 
the different shapes of their power spectra. 

As stressed by Scott & White (1999), the amplitude of 
the clustering signal drops rapidly - in comparison to the 
one arising from primary anisotropies - with decreasing fre- 
quency. In the nearest Planck/HFI channel (centered at 217 
GHz) its contamination on the CMB anisotropy maps should 
already be negligible. 
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